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Abstract The partition function of the two-dimensional Ising model on a square lattice with nearest-neighbour 
interactions and periodic boundary conditions is investigated. Kaufman [Phys. Rev. 76, 1232-1243 (1949)] 
gave a solution for this function consisting of four summands. The summands are rewritten as functions of 
a low-temperature expansion variable, resulting in polynomials with integer coefficients. Considering these 
polynomials for system sizes 2™ x 2" {in,n G N), a variety of recursion relations in in,n are found. The re- 
cursions reveal a rich structure of the partition function and can be employed to render the computer algebra 
calculation of the microcanonical partition function more efficient. 
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1 Introduction 

Every scientific community has its favourite model system, and the role which is played in genetics by 
drosophila melanogaster can be attributed to the Ising model in the statistical physics of phase transitions. 
It is arguably the best studied model in the field, and advances in the understanding of the Ising model have 
significantly advantaged progress in the study of phase transitions in general. 

Proposed by W. Lenz in 1920, the one-dimensional version of this lattice model with nearest-neighbour 
interactions was solved by his student E. Ising |3 1 in 1925. Ising showed the absence of a phase transition in 
the one-dimensional case, and this result, enforced by king's claim that likewise no phase transition should 
take place in higher dimensions, discouraged further studies for a while. It was Peierls' proof |7 2| of the 
existence of a phase transition in the two-dimensional Ising model on a square lattice which rekindled the 
interest, culminating in Onsager's exact solution |6l for the free energy density, published in 1944 (see for 
an account of the early history of the Ising model). 

The above cited articles all focus on properties of the Ising model in the thermodynamic limit, in which 
the number of lattice sites goes to infinity. In 1949, B. Kaufman published an extended version of Onsager's 
solution, reporting an exact expression also for the canonical partition function of the two-dimensional Ising 
model on a finite square lattice with periodic boundary conditions |4|. Almost 50 years later, P. D. Beale 1 1 1 
employed a computer algebra program to calculate from Kaufman's result, for modest lattice sizes of up to 
32 X 32 sites, the exact microcanonical partition function for the energy, a quantity particularly useful when 
testing the performance of computer simulation algorithms. The microcanonical partition function, taking on 
values from the non-negative integers No, is given by the number of configurations compatible with a certain 
value of the energy. From Beale's algorithm, these integer numbers are computed from sums and products of 
irrational numbers. This does not appear to be the straightforward way to calculate an integer number, and it 
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was this observation which made me start investigating the structure of the finite system partition function of 
the two-dimensional Ising model. 

When working on this problem, I came across beautiful numbers and rich structures, mainly in the form 
of recursion relations between polynomials with integer coefficients. Although the solution obtained is not 
complete (in a sense which will become clear later), I would like to share at least the problem, its beauty, 
and a partial solution with the scientific community. Apart from the esthetic appeal, the study can have the 
practical application of improving the performance of Beale's algorithm. 



2 The two-dimensional Ising model 

Consider a finite two-dimensional square lattice 

^M,yv = {l,...,M}x{l,...,7V}, (1) 

consisting of MN lattice sites. Periodic boundary conditions are assumed, and to each lattice site ; G J^m.n a 
classical spin variable a, G { — 1, +1} is associated. Then the two-dimensional Ising model (without external 
magnetic field) is defined by the Hamiltonian function 

^(a) = -y £ a,a„ (2) 
mapping each configuration a = from the configuration space 

rMN = {'\.+\Y'''' (3) 

onto J times the even integer numbers 

{0,2,4..., ±2MN - 6, ±2MN - 4, ±2MN} . (4) 

The positive constant / determines the coupling strength and the angle brackets ( • , • ) in (|21l denote a summation 
over all pairs of nearest neighbours on the lattice J^m,n (each pair counted once). 



3 Kaufman's expression for the canonical partition function 

The canonical partition function is given by 

^M^(K)= £ e-^^'-f-), (5) 

where K = 2/j3 and j3 = ^ is the inverse of temperature T times Boltzmann's constant kg. The ^T-dependence 

of will be suppressed in the following to lighten the notation. Kaufman |4l showed that the canonical 
partition function of the two-dimensional Ising model can be written in the form 

where 

= {IsinhKf/' n 2cosh (^Y2k+i) , (7a) 

k=0 \ ^ / 

J-^^ = (2sinh7^)^^/2 n 2sinh f f y^^+i) , (7b) 

k=0 V ^ / 

3fll = (2sinh^)'^^/2 n 2cosh (^72^ , (7c) 

k=o \ ^ / 

iT^) = (2sinh^)^^/2 n 2sinh (^y^,) . (7d) 

k=0 V ^ / 
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For < < 27V the quantity J]^ is defined as the positive root of 



sinh/: V N 

whereas for = we have 



cosh(}^) = ^"^"f 5 - cos f ^ ) , (8a) 



e 
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= e^tanh(7i:/2). (8b) 



The invariance of upon interchange of the parameters M and A', ahhough evident from the definition of 
the Ising model, is not quite obvious in the above expressions. 

As a consequence of the equidistant energy values the canonical partition function can be expressed as 
a low-temperature series of the form 

MN 

^MM=^-''''j:8kx"', (9) 

k=0 

where 

x = e-^, (10) 

is the low-temperature expansion variable. This is exactly what is done by the computer algebra program at 
the core of Beale's work 1 1 1: For given M and A', the solution specified by equations (|6|l-(|8j is cast into the 
form (|9). From this result the microcanonical partition function gj^, specifying the number of configurations 
with energy —2J{MN — 2k), can be read off. The gj^ are non-negative integers, but their computation from 
equations (|6}-(|8} involves sums and products of irrational numbers, stemming essentially from the cosine in 



4 Low-temperature expansions of the 

Since low-temperature expansions show the appealing property of having integer expansion coefficients, we 
continue working in this representation. One can 
(|6j in the form of a low-temperature expansion, 



continue working in this representation. One can similarly express the four summands (i = 1,2,3,4) in 



with expansion coefficients g^^ . From Kaufman's solution (|6|l-(|8} or from Beale's computer algebra program 
in yj one obtains for example 

(12a) 
(12b) 
(12c) 



(12d) 



-^-4 

X ^2,2 




)^ 




^16 ^(1) 
^4 4 


=(1+.- 


)'(l-4x2 + 


22/-4/+jc^)^ 


^36 avW 




)'\\-%^^ 


- 30/ - 8/ +.x**)'(l +/ + 12/ +/ +/)4, 


^ °^8,8 


=(i+x2; 


)'''(l+20/ 


+ 32/ + 150/ + 32jc1° + 20jc'2 + jc^^)'* 




x(l- 


8^2+28/- 


56/ + 326/ - 56x"' + 28jc'2 ^ 3^14 _^^i6^ 2 



Remarkably, not only the g]^ from the series (|9j, but already the in ( fTTT i are integer numbers (not nec- 
essarily non-negative, however). So the S^'^-^ are essentially polynomials in x with integer coefficients. The 
complexity of these polynomials can be reduced significantly by the transformations of variables 

(\-xW 4(l+v)2 
y=\-^) and ^= ' (13) 
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yielding for example 

^2,2 =(4^)' 2, (14a) 

^i,P=(4}')'z'(8-z)', (14b) 

2f^^={^yfz\\2-zf{^-z)\ (14c) 

J^^(i) = (4);) (32 - 16z + z') ' (8 - 8z + z') ' . (14d) 

In the following we will establish recursion relations between the canonical partition functions and 

2^214 N '^^ "^M 2N' respectively. To this purpose, and for further simplification, we restrict ourselves to the Ising 
model on lattices ^2"" ,2" where m,« G N. Then it is convenient to introduce the notation 

^m,n ~ (4j) ^2^,2'^ ^ (15) 

yielding for example 

Zia' =z, (16a) 

Z^=z2(8-z)^ (16b) 

z\^l =z^ {32 - 16z + z^f {8- 8z + z^)\ (16c) 



4 4 (512 - 1024z + 320z2 - 32z^ +z'*)^ (8 - 32z + 40z2 - 16z^ 



X (8 - 96z + 72z2 - 16z^ + z'^) (32 - 128z + 80z^ - 16z V z^) \ (16d) 

(compare with equations ( I14a> . ( I14b> . and ( I14d» . The rest of this article is devoted to the study of relations 
among the Zm,n- Instead of having tried to obtain these relations from the general equations (|6}-(|8}, the 
Zmji were calculated by means of a computer algebra program for several values of m and n from which 



the relations were "read off". It is for this reason that no derivations are given in sections |5H7] The validity 
of the relations has been verified for values of m,« up to 7, corresponding to linear lattice sizes up to 128 
sites. The present article is supplemented by a MATHEMATICA notebook available at http://www.phy.uni- ^ 
bayreuth.de/ btp308/Zmn.tar.gz containing a data set and some routines which allow to construct and handle 

the Zm'n, given as polynomials in z. 



5 Symmetry relations for the Z„,.„ 

For any given m,n G N, the symmetry relations 

7(1) _ 7(1) 7(4) _ 7(4) . 7(2) _ 7(3) .,7. 

hold. In combination, these relations imply the invariance 

3f2'",2" = ^2",2'" (18) 



of the canonical partition function. 



5 



6 Relations between the zlj^„ 

The Zm,n are polynomials in z of degree 2™+""^. As may already be suggested by ( fT6i . these polynomials 
typically factorize into several polynomials of lower degree. Many of these factors appear repeatedly in the 

various Zm'„, and, for all m,n 2, the relations 

7(3) _ 7(1) y(3) _y(2)rT7(l) .jq.. 

ZW=ZP) ,Z(^) ,=z|4)nz*'l, (19c) 

m.n—l m.n—l l.m 1 1 m.A:' \ ^ / 

can be read off. To compute the various Z^]„ («' = 2, 3 , 4) from these relations, the z{'^ = 2, 3 , 4) are needed 
as starting values, for which we find 

ZS=^, (20a) 

Z?i=biUal (20b) 

The fl„ are given by the recursion relation 

a\ = v^, (21a) 

for n ^ 2, (21b) 



yielding a sequence whose first few terms read 

ai = ^/z, (22a) 

a2 = ~2 + z, (22b) 

a3=2-4z + z^ (22c) 

c^4 = 2-16^ + 20^2-8^^ + z^ (22d) 

fl5 = 2 - 64z + 336z2 - 672z^ + 660z^ - 352z^ + 104z^ - 16zV z**. (22e) 
The bn are determined by 

bi=-4 + z, (23a) 

Z?„ = (^„-i+2)2-8j3„_i forn^2, (23b) 

where 

Pi=z, (24a) 

j8„ = |^j3„_i - 2 + 4«-izn j3,^ for n ^ 2, (24b) 

and the resulting sequence starts with 

bi = -4 + z, (25a) 

Z72=4-12z + z^ (25b) 

/73=4-176z+148z^-24z^+z^ (25c) 

b4=4- 2752z + 29520z2 - 52704z^ + 30356z'^ - 7456z^ + 872z*' - 48z^ + z^ . (25d) 
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The resuhs of this section effectively reduce the calculation of the z'^^n for any i,m,n to the computation of 
theZ^') 



7 Recursion relations for the z/J i 

For an efficient computation it would likewise be desirable to find recursion relations for the zjnli for different 
m^n. For = 1,2 (or, equivalently, n = 1,2), Z/fi^/i can be written as 

Z['l=al and 4m=<^1 (26) 

with the a„ as given in (121 > . The c„ are, again, given by recursion relations, 

ci=-2 + z, (27a) 
C2=z(z-8), (27b) 

c„ = (c„_i+2)2-l(c„_i-c2_2+4)2 forn^3, (27c) 

yielding a sequence that starts with 

ci=-2 + z, (28a) 

C2 = z(-8+z), (28b) 

C3 = 4-32z + 60z^-l6z^ + z\ (28c) 

C4 = 4 - 128z + 1264z2 - 4160z^ + 46282"* - 1984z^ + 376z^ - 32z^ + z^ . (28d) 

(1)' 



The next sequence of interest would be {Z3 j starting as 



4]l={2-4z + z^)\ (29a) 



(4 - 32z + 60z2 - 1 6^3 + z"*) ^ (29b) 

z[]l =z^(32-16z + z2)^(8-8z + z2)\ (29c) 
4a = (4 - 352z + 3424z2 - 6976z^ + 5436z'* - 2000z^ + 368z^ - 32z^ +z^f 

X (4 - 160z + 1376z^ - 3648z^ + 3900z'* - 1840z^ + 368z^ - 32z^ + z^)^ (29d) 



-"3,2 
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^3.5 =(4 - 512z + 23936z^ - 504128z^ +4799728z'* - 20395520z^ 

+ 44112832z^ -55088640z^ +42942996z^ - 21874432z'^ +7464768z 

- 1722784z^' +267640z^^-27392z^^ + 1760z''*-64z'^+z'^)^ 
X (4 - 1536z+ 149888z^ - 1832128zV9714928z'* - 28345856z^ 

+ 50377664z^ - 57468672z^ +43287060z^ - 21857536zV7456064z^° 

- 1722208z^' +267640z'^ - 27392z^^ + 1760z''* - 64z'^ +z'^)^ (29e) 

(for more terms of this sequence see the accompanying MATHEMATICA notebook at |http://www.phy.uni-| 
bayreuth.de/ btp308/Zmn.tar.gz). I was not able to obtain recursion relations (in n) for this sequence, neither 

for the other Z^j, with in ^ 4. Similarly, recursions for zi]),, i.e., for the Ising model on square lattices of 
linear size 2", would be desirable. 
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8 Summary 

We investigated Kaufman's solution |4| of the partition function of the two-dimensional Ising model on an 
M X N square lattice with nearest-neighbour interactions and periodic boundary conditions. The four sum- 

mands 5^''^ {i = 1,2,3,4) contributing to this solution were considered separately. Rewriting these sum- 
mands as functions of the low-temperature expansion variable x = exp(— 2/j3), polynomials with integer 
coefficients were obtained. The polynomials can be simplified significantly by the change of variables defined 
in il3\ . Restricting the analysis to lattices of sizes 2'" x 2" {m,n G N), a variety of recursion relations in m,n 

were found for the resulting simplified polynomials, and hence for the ^j^i^- These recursion relations re- 
duce the calculation of the Zm?„ for any ;,ot,« to the computation of Zmj,. For m= 1,2, the z/,,], can again be 

computed efficiently from recursion relations for arbitrary n. The search for similar recursions for Zm j, when 
m ^ 3 was unsuccessful and remains an open problem. 

The results obtained are of interest for different reasons. First, the recursion relations can be employed to 
render more efficient the (computer algebra) calculation of the microcanonical partition function in the spirit 
of Q. With the results obtained so far, the gain in efficiency is still modest; if, however, recursion relations 
like the ones indicated as open problems in sectionQcould be found, the efficiency of the computation would 
increase drastically. 

Second, the results reported in this article reveal a fascinatingly rich structure of the partition function. 
In my opinion it is in particular the reformulation of the ^^''^ as polynomials in z with integer coefficients 
which contributes significantly to the esthetic appeal of the problem discussed in this article. My hope would 
be that the results presented attract the interest of someone capable of making further progress in the direction 
pointed out. 
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